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Dynamical Stability of Cusps and Cores: 
Implications to Centers of Galaxies and Clusters of Galaxies 

Tadashi Nakajima 1 *) **) 

1 National Astronomical Observatory of Japan, Osawa 2-21-1, Mitaka, 181-8588, 

Japan 

We study the stability of a spherically symmetric density profile. First, we examine 
the case for ordinary gas with a well-defined equation of state (EOS). In the Newtonian 
regime, the analysis is based on the equation of hydrostatic equilibrium, and in the general 
relativistic regime, it is based on the Tolman-Oppenheimer-Volkoff equation. We find that 
a stable system must have a flat top, or dp/dr — at r = 0. For a power-law equation 
of state, p oc p 13 , (3 > is required in the Newtonian regime, while f3 > 1 is necessary 
in the relativistic case. At least, some contribution from degeneracy pressure of fermions 
(1 < j3 < 5/3) is necessary to support a stable structure in the case of relativistic gas. In these 
cases, stability conditions are justified thermodynamically as well as dynamically. Second, 
we analyze the case of a collisionless system with a power-law profile given by p oc r~ a , 
in the Newtonian regime using the Jeans equation. The Jeans equation is mathematically 
identical to the equation of hydrostatic equilibrium as long as p.j — pa 2 is regarded as the 
pressure counterpart for ordinary gas, where a is the local velocity dispersion. Since a self- 
gravitating collisionless system is at best metastable, if it is treated statistical mechanically, 
the dynamical equilibrium described by the Jeans equation will have a shorter life time 
than the thermodynamically stable equilibrium for ordinary gas. Here we take a heuristic 
approach, using the pressure counterpart, pj. From the requirement that pj > or the 
normalcy of the formal EOS, pj = pj(p) obtained by eliminating r, we find the following 
heuristic stability criterion. For < a < 1, pj < 0, the density profile is unstable. We 
interpret this profile to be a transient one evolving toward a flat-top profile, which seems to 
be the final product with absolute stability. For 1 < a < 3, the profile is dynamically stable. 
It is an interesting coincidence that dark matter halo profiles obtained by cosmological N- 
body simulations, fall into this category. For 3 > a, pj < 0. We suspect that such a 
steep profile is in the process of collapsing toward r = 0. Finally, our stability criteria are 
compared with the density profiles obtained for galaxies and clusters of galaxies. 



§1. Introduction 

For ordinary gas, with a well-defined equation of state (EOS), the equilibrium 
structure of a spherically symmetric self-gravitating system such as a star, can be 
obtained by integrating the equation of hydrostatic equilibrium (HSE) in the New- 
tonian regime and the Tolman-Oppenheimer-VolkofF (TOV) equation 5 ) in the rela- 
tivistic regime. In this process, both inner and outer boundary conditions must be 
specified. At the inner boundary, both the pressure, p, and density, p, need to be 
positive and finite. By integrating the equation of HSE or TOV equation, we show 
that these requirements on p and p lead to the requirements on the flat-top nature, 
dp/dr = at r = 0, of the density profile and on the power-law index of the EOS, 
(3 > 1, for p oc pP. A physical interpretation of the requirement on (5 is that at least 
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some contribution of degeneracy pressure of fermions as a repulsive force is necessary 
to support a dynamically stable structure. In the real universe, neutron stars, white 
dwarfs, and cool brown dwarfs, are supported by degeneracy pressure of neutrons or 
electrons. 

Now we turn to the analysis of a spherically symmetric system composed of 
collisionless particles in the Newtonian gravity. ^ There is mathematically a formal 
analogy between the equation of HSE and the isotropic Jeans equation, if we define 
the pressure counterpart pj = pa 2 , where a is the local velocity dispersion. The 
Jeans equation describes the situation in which the radial gradient in the product of 
the density and velocity variance is balanced with the gravity induced by the encircled 
mass at a given r. As is well known, statistical mechanics of a self-gravitating system 
is not a well-defined subject due to the long-range nature of the Newtonian gravity. 6 ) 
The equilibrium structure described by the Jeans equation is at best metastable 
thermodynamically. On the other hand, the dynamical equilibrium is realized within 
a much shorter time scale than in a thermodynamical relaxation time, and the Jeans 
equation is expected to describe the stability on a time scale much shorter than the 
thermodynamical stability. 

To proceed with the analysis of the stability of a density structure described by 
the Jeans equation, we make a couple of assumptions regarding the analogy between 
the equation of HSE and the Jeans equation. First, the inner boundary condition 
imposed on the EOS for the ordinary gas, is the same as that for the formal EOS 
introduced for pj for a collisionless system. Second, there is an analogy between the 
thermodynamical normalcy of the gas, which appears in the behavior of the EOS, 
and that in the formal EOS for the collisionless particles. For instance, we regard a 
negative pj as an indicator of strong instability of the apparent equilibrium described 
by the Jeans equation. 

We then apply this approach to a structure described by a power-law density 
profile given by p oc r~ a . At a finite r, elimination of r in the Jeans equation gives an 
"EOS" or the relation between pj and p. If our criterion turns out to be correct, it 
becomes possible to describe the strength of "met ast ability" , or "instability" , which 
depends on the value of the power-law index, a. 

We admit that these assumptions are not fully backed up by a sound physical 
principle. Our justification for now, is more historical than scientific. In the past, 
some heuristic approaches proceeded to the development based on rigorous physical 
principles. What we hope here is that if we find some sound physical basis of this 
approach in the future, our stability criterion will become more reliable, since the 
results of mass distribution measurements of astronomical objects and cosmological 
simulations, are often expressed by this power-law form. 

The paper is organized as follows. We first elaborate the derivation of the 
stability criteria for ordinary gas in §2. Our main discussion on the heuristic stability 
analysis for collisionless particles is given in §3. The applicability of our stability 
criterion to centers of galaxies and clusters of galaxies, and density profiles obtained 
by N-body simulations, is discussed in §4. 
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§2. A Necessary Stability Condition for a Self-Gravitating Gaseous 

System 

2.1. Non-Relativistic Regime 

The basis equations are the equation of hydrostatic equilibrium (HSE), 

■m _ , (2 .i) 

combined with the encircled mass relation, 

alM{r) = Airr 2 p(r)dr. (2-2) 

In the case of a star without nuclear energy generation, these equations are simul- 
taneously integrated, starting from either the inner boundary or the outer boundary, 
for a given equation of state (EOS), p = f(p). For a dynamically stable structure, 
there is a necessary condition imposed at the inner boundary at r = 0, that p is 
positive and finite at r = 0. 

Now we consider the case that the EOS is given by a power-law as, 

p = AfJ 3 . (2-3) 

As long as (3 is finite, p, is also finite at r = 0, which gives the boundary condition 
for M(r) as M(0) = 0. 

Ea (|2-lj) for a small r can be expressed as, 

^ = -4GW0) 2 r. (2-4) 
ar 

Then, dp/dr — > as r — > and dp/dr < for a finite r. From the same 
argument, we obtain the similar behavior of p for positive A and (5. 

Therefore, the stability near the origin requires the following two conditions, 

(1) The density structure must have a flat top. 

(2) For the EOS, p and (5 must be positive. This means that p increases as p increases, 
indicating thermodynamical normalcy. 

2.2. General Relativistic Stability and Equation of State 

At extreme high densities near the origin of symmetry of mass distribution, 
we have to take into account at least two physical phenomena outside the scope 
of the classical Newtonian mechanics, the quantum statistical degeneracy, and the 
general relativistic treatment of stability. Even if dark matter particles are only 
weakly interacting in the sense of particle physics in vacuum, indistinguishability 
of identical particles affects statistical dynamics, when the mean particle spacing 
becomes smaller than the thermal de Broglie wavelength. 

Actually this is not exceptionally remote possibility. At the core of a massive 
cluster of galaxies, the mass density reaches 10 -24 g cm -3 . If a significant fraction 
of dark matter is composed of eV-mass massive neutrinos, they are expected to be 
degenerate on the scale of 100 kpc. 2 ) 
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Although a collisionless Newtonian system composed purely of cold dark matter 
may be thermodynamically metastable, the effect of repulsion by degeneracy pressure 
of neutrinos may be to make the mixture of neutrinos and cold dark matter particles 
into a thermodynamically normal system. 

Here we assume that the equation of state (EOS) of the dark matter gas will 
be well defined in the local inertial frame near the origin. The EOS in terms of the 
pressure law, p = p(p), governs the behavior of the energy momentum tensor, and 
thus the stability. Again we consider the EOS expressed as a power law as, 

p = Ap^, (2-5) 

where A does not depend on p, but may depend on temperature, T, and particle 
properties such as the particle mass and statistical weight. We examine the range of 
the power law index, (3, for bosons and fermions, under different physical conditions. 
First we list up the extreme cases and then give the ranges of (3 for bosons and 
fermions separately. 

For nonrelativistic and non-degenerate gas (Maxwell-Boltzmann ideal gas), 



13 = 1. (2-6) 

For extremely relativistic bosons, 

= 0, (2-7) 

since they can be regarded as a black body for which particle numbers are not 
conserved. 

For nonrelativistic partially degenerate bosons, most particles are in the ground 
state and remaining thermally excited particles contribute to the pressure, p. p 
depends on the number of thermal particles determined by temperature, T, but not 
by p or the total number of the particles. 4 ) 



(3 = 0. (2-8) 

Including the cases with intermediate physical conditions, the range of (3 for 
bosons, is 



P = ~ 1. (2-9) 

In the case of fermions, the range of (3 for fully degenerate ideal gas is between 
4/3 (extreme relativistic) and 5/3 (nonrelativistic). 4 ) Including the cases with inter- 
mediate degeneracy, the range of (3 is 



/3 = l~5/3. (2-10) 

This problem setting under general relativity may appear formidable. However 
it has been dealt with by the classical paper by Oppenheimer and Volkoff (OV39) 
on the neutron stellar core. 5 ) They first formulated the problem based on Tolman's 
coordinate system for spherical symmetry 7 ) for a general equation of state and then 
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applied the particular parametric equation of state for fully degenerate fermions 
covering nonrelativistic to extreme relativistic particles 4 )' 8 ) to solve the neutron star 
problem, assuming neutrons as ideal fermi gas. The general relativistic counterpart 
of the equation of the hydrostatic equilibrium is the TOV equation after Tolman, 
Oppenheimer and Volkoff. For the generalized radial coordinate r in the Tolman's 
coordinate system and corresponding pressure gradient dp/dr, OV39 discuss the local 
stability requirement in the footnote. The result is that the local stability is satisfied 
for /3 > 1. 

Here we briefly summarize and elaborate the argument by OV39. Using the 
Tolman's coordinate system and the dimensionless encircled gravitational mass (en- 
ergy), -u(r), Einstein's equations are reduced to two equations, 

^ = 47re(p)r 2 , (2-11) 
dr 

and 

Art m -I- £\ r\ 1 „ 

-[Anpr 3 +u], (2-12) 



dr r(r — 2u) 

where e(p) is the energy density and p is the pressure in the proper coordinate 
of the matter. OV39 require that at r = 0, u and p must have uq = and a positive 
finite value, po respectively. If e(p) is given as a power law of p as, 

e(p) = C t p s , (2-13) 
near the origin, ea (j2T2|l reduces to 

dp p + C t p s 



dr 2r 

For po to be positive and finite, integration of ea ()2T4|) shows that s < 1. 
Now we examine the relation between the power- law EOS, p = Ap@, and e 



(2-14) 



pc 2 + e K , (2-15) 



where ex is the kinetic energy density. 

In the nonrelativistic case, pc 2 dominates e and 



c» 4) . (2-16) 



.A 

So one can see that s < 1 is equivalent of (3 > 1. 

In the ultrarelativistic limit in which gas is treated as radiation, 

P=\e, (2-17) 

or s = 1. Therefore there is no stable solution for ultrarelativistic case. One 
may worry that the Chandrasekhar limit corresponds to the EOS of ultrarelativistic 
degenerate electrons in white dwarfs. In the case of white dwarfs, the gravity source 
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is hadrons, which are nonrelativistic and gravity is Newtonian. On the other hand, 
the TOV equation is dealing with "self-gravity" of the matter. Even at the center 
of a neutron star, neutrons are only moderately relativistic. 

We need to examine more carefully for the moderately relativistic case. In 
the nonrelativistic limit, for ideal gas regardless of bosons, fermions or degree of 
degeneracy, 

V=\zk- (2-18) 
and from the relativistic limit, eq (|2-17|) . 

P = \e K <\e. (2-19) 

So it is apparent that the second term of ea (|2-15|) . ek alone never prevents the 
structure from collapsing at the origin. So the stability at the origin is fully deter- 
mined by the behavior of p as a function of p. Therefore (3 > 1 is always required. 
For j3 > 1, dp/dr = at r = and the integration of the TOV equations gives a 
flat-top density profile again at a small r. To summarize, the stability condition near 
the origin requires, 

(1) The density structure must have a flat top. 

(2) For the power-law EOS, p > and j3 > 1 are required. Since the system 
is unstable for (3 = 1, classical thermal pressure can never support a relativistic 
structure. To form a stable relativistic structure, some contribution from degeneracy 
pressure by fermions, or the repulsive force due to Pauli's exclusion principle, is 
required. 

§3. Stability of a Power-Law Density Profile for Collisionless Particles 

in the Newtonian Regime 

3.1. A Heuristic Treatment of the Jeans Equation 

Now we turn to our heuristic analysis of a spherically symmetric collisionless 
system in the Newtonian regime. In the following, we use the Jeans equation, trying 
to find some analogy with the equation of hydrostatic equilibrium for ordinary gas. 
We formally introduce the pressure counterpart and the counterpart for the EOS. 
As is well known, a self-gravitating system formally treated by the micro-canonical 
ensemble of statistical mechanics, is thermodynamically metastable, when the con- 
tributions of the kinetic energy and gravitational potential are comparable. 3 ^' 6 ) This 
is exactly the situation that the Jeans equation is dealing with. So first, we accept 
this thermodynamical anomaly, but expect some type of dynamical stability, when 
the Jeans equation is valid. Below we show that depending the steepness of the 
power law, the pressure counterpart becomes negative and the formal EOS indicates 
strong anomaly. In this case, we question the applicability of the Jean equation itself 
and the corresponding density profile is regarded as "strongly metastable" or even 
"dynamically unstable", compared to a system with normal EOS under the Jeans 
equation. We also caution that our "pressure" or EOS approach is not same as the 
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treatment of the statistical mechanics either. To define a pressure, in the sense of 
standard statistical mechanics, we must define a boundary to define a volume. How- 
ever, this procedure is inconsistent with the long-range and unshielded nature of the 
gravity, and therefore "pressure" of this sense was never defined by the approach of 
statistical mechanics. 

3.2. A Power-Law Density Profile 

The results of mass distribution measurements of astronomical objects, such as 
those of galaxies and clusters of galaxies, are often expressed as density laws in the 
form of 



In the following, we examine how much information is extracted from this phe- 
nomenological profile. We begin with some certain facts, and then proceed to our 
heuristic analysis of dynamical stability. 

3.3. Encircled Mass and Velocity Dispersion 

The most simple physical quantity for a given density law, ea (j3-l|) . is the encir- 
cled mass, M(r). 



So M(r) — ► for r — > 0, only if a < 3. However M(0) = is only a necessary 
condition for a real astronomical object, since even a black hole has a finite extent. 
If one finds a profile with a > 3 in a certain range of r, it indicates the presence of 
an extra-mass source inside, indescribable with the power-law alone. 

The second quantity to consider is the velocity dispersion, c(r), which is often 
an observable. 



Then, <r(r) — ► oo as r —* for a > 2. In reality, the relativistic requirement, 
a < c, sets the limit on the minimum value of r m i n , which gives an order estimation 
of the Schwarzschild radius, 2GM/c 2 . So if a cusp shows a rise of density towards the 
center with a > 2, this is an indication of the presence of a compact source inside. 
It may appear that there is no anomaly for a < 2. However, a simply indicates the 
behavior of test particles and it is not necessarily related to the stability. 

3.4. Jeans Equation 

For an isolated spherically symmetric system with collisionless particles, such 
as that of cold dark matter particles or an idealized stellar system with equal mass 
stars, the dynamical stability condition is expressed by the isotropic Jeans equation 1 ) 








(3-3) 
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as, 

d[n(r)a 2 (r)} [n(r)GM(r)] 
dr r 2 

with 



(3-4) 



M(r) = / A-Kmn(r)r 2 dr, (3 - 5) 
Jo 

where n{r) is the particle number density, m is the particle mass. 
Suppose we set 

p(r) = mn(r) (3-6) 

and 

pj = mn(r)a(r) 2 , (3-7) 

where we note that a{r) is the radial component of the total velocity dispersion, 
°"toi = \/3cr, then the Jeans equation and the encircled mass can be rewritten as 



dp j GM(r 

— = -p(r) ; 

dr r 

and 



— = -P\ r > ^2 » ( 3 ' 8 ) 



M(r) = / 4irp(r)r 2 dr. (3-9) 



o 

Here we caution that an observable velocity dispersion is often the component 
in the line of sight, a, instead of otot- 

Therefore, the Jeans equation for collisionless particles and the equation of hy- 
drostatic equilibrium are identical, as long as we regard 

pj = mn(r)a(r) 2 = p(r)a(r) 2 , (3-10) 

of a collisionless system as a counterpart for the pressure, p of ordinary gas. This 
analogy is at least energetically reasonable, since the pressure is two thirds of the 
kinetic energy density for nonrelativistic ideal gas. However from the point of view 
of particle momenta or force, there is a clear distinction. Gas pressure is the total 
momenta given to a boundary wall per unit time, while no boundary can be defined 
for the self-gravitating system. So we should regard this analogy as mathematical for 
now. Now we examine the behavior of the power-law profile by the Jeans equation 
and pj in detail. 

If we integrate ea (|3-4j) . for the power-law profile, ea (|3-lj) . 

Pj {r) = A p r 2 ^+ Pj (0) (3-11) 

and 
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where A p is a constant coefficient and p(0) is a constant of integration. When 
p(r) is given as ea (|3-l|) . a(r) is also another power law as ea ()3-3|) . then pj must also 
be a power law. Therefore pj(0) = 0. So 

pj(r) = A p r 2 ^ a \ (3-13) 
where r 2 ( 1-a ) > for a finite r, and the sign of pj is determined by A p as, 



A p < for a > 3, (3-14) 
> for 1< a < 3, (3-15) 
< for a < 1. (3-16) 

Both p(r) and pj(r) are power-laws of r and by eliminating r, we obtain a formal 
EOS as, 

pj = A p p\ (3-17) 

where, 

7 = 2(l-i). (318) 

Now we examine the usage of the pressure counterpart pj and the EOS coun- 
terpart, ea (|3-17|) . 

3.4.1. Thermo dynamical Instability at the Origin 

From the result of §2, thermodynamical stability at r = requires, pj to be 
positive finite and the density profile to be a flat top. The former is satisfied only 
for a = 1, but the density profile, p(r) oc r" 1 , does not have a flat top. So there is 
no power-law profile that satisfies "thermodynamical stability" . If a flat-top at the 
origin is combined with a smoothly dying density profile near the outer boundary, 
such as that for a star, there must be an inflection point of the density profile, 
somewhere at an intermediate r. On the other hand, a single-power-law profile can 
never have an inflection point and the situation is understandable. 

3.4.2. The behavior of the Pressure and EOS Counterparts at a Finite r 
At a finite r, if we require the positivity of pj, a must be in the range, 



1< a < 3. (3-19) 

If we further require that the EOS is normal, in other words, pj increases as p 
increases, 7 > 0. It turns out that this requirement also gives, 1 < a < 3, equivalent 
to the positivity condition for pj. Since a system described by the Jeans equation is 
at most met ast able or in dynamical equilibrium with a finite life time, we consider 
that pj < 0, indicates even shorter life time of this stage. 
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3.4.3. Some Conjecture About Physics of Pressure Counterpart 

For a > 3, from ea (|3-2p . M(r) — > oo as r — > 0. This is an unphysical situation 
and a single power-law profile cannot describe a stable structure. In that sense, 
pj < 0, may indicate some collapsing stage to the origin. Or pj < can be an 
indication that the validity of the static Jeans equation itself is violated, indicative 
of dynamical instability. 

The range of a = 1 ~ 3, for positive pj, has an interesting coincidence with 
the halo density profiles obtained by cosmological N-body simulations. 10 ^ 13 - 1 For 
instance, the original NFW profile behaves as a — > 1 at small r, while a — > 3, as 
r — > oo. From the encircled mass argument based on ea (|3-2j) . a finite total mass at 
r — > oo requires a > 3 at a large r, and the behavior at a large r, does not seem to 
carry much dynamical information, a ~ 3, also indicates that the outer region of the 
NFW profile corresponds to pj ~ 0, indicative of invalidity of the Jeans equation. 
So our main interest is in the stability of cusps at a small r, obtained by N-body 
simulations. Including other simulations, the power-law indices of inner cusps are in 
the range of a = 1 ~ 1.5. This range of a should show a dynamical equilibrium with 
a long life. 

The most interesting region of a is a = ~ 1. As we have seen in §2, a 
thermodynamically stable profile must have a flat-top. The density profile of a star, 
for, instance, has a flat core, indescribable with a single power law. If the flat core 
profile is the thermodynamically stable final stage, it appears reasonable to regard 
this range of a to indicate the transition from the dynamically stable cusp with 
a = 1 ~ 1.5 to a thermodynamically stable flat-top core. 

The coarse grained view adopted by N-body simulations will give a proper time 
scale in the free-fall stage, before dynamical friction sets in. However once dynamical 
friction becomes significant at a high density region, the time scale for dynamical 
evolution behaves as N/ log iV, 6 ) where N is the total number of collisionless particles 
and the particles number or particle mass must be known to evaluate the evolution 
time scale. So the actual duration of this possible transient stage may not be esti- 
mated until we know the mass of a dark matter particle, even if our conjecture is 
correct. The summary of our stability arguments is given in Table [I] 



Table I. Thermodynamical and Dynamical Stability of Density Profiles Based on Our Criteria 



Profile Type 


a 


Thermodynamically 


Dynamically 


PJ 


comment 


Flat Top 


NA 


Stable 


Stable 


Positive 


Final Stage? 


Power Law 


< a < 1 


Unstable 


Unstable 


Negative 


Transient? 




1< a < 3 


Unstable 


Stable 


Positive 


NFW in this range 




Q > 3 


Unstable 


Unstable 


Negative 


Collapsing to origin? 
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§4. Current Observational Status of the Core/Cusp Problem 

4.1. Limitations in Observations and Applicability of Our Stability Criteria 

Before comparing the observational status with our heuristic stability criteria, 
we must carefully examine the limitation in observations to avoid over interpretation. 

4.1.1. Indistinguishability between a Flat-Top Profile and an Unstable Shallow Cusp 
The spatial resolution of observations is always limited and even if the true 

density profile is not a single power law, the observational results are often fit by a 
single power law for a limited range of r. For instance, even if the true profile is a 
flat-top type with an inflection point, one can fit its observed profile given for a finite 
range of r with observational errors, using a single power-law, and possibly find a 
shallow cusp with < a < 1. The real observational situation seems to be that we 
are not able to distinguish a flat-top profile from a shallow cusp for external galaxies 
or clusters of galaxies. 

4.1.2. Dark Matter, Baryons and Stars 

The density profile of a spiral galaxy is usually obtained by the measurements 
of a rotation curve of the interstellar medium. What rotation curve gives is the 
encircled mass inside the given radius. On the average over the entire universe, the 
mass ratio between the dark matter and baryons is about 5:1, but we do not know 
of this true ratio, especially at the inner region of a galaxy. The actual presence 
of stars is a manifestation that the dissipative collapse and star formation actually 
happened in the past. Therefore there is no guarantee that the dark matter and 
baryons or stars follow the same density profile. 

One possible relief to this complicated situation is that due to the identity of 
the mathematical form between the equation of hydrostatic equilibrium and the 
Jeans equation, a similar stability argument may be possible even for the mixture 
of ordinary gas and collisionless particles, as long as we know the total density 
profile, measured by the rotation curve. However this is based on a very optimistic 
assumption that the dynamical time scales for dark matter, baryons, and stars, 
are the same, which is not necessarily true. Unless the entire system is already in 
dynamical equilibrium, this treatment will be invalid. It is probably wrong, if we 
analyze the stability of the system with a mixture of various types of particles, using 
our criteria. 

This situation seems severe for individual galaxies, while for a cluster of galaxies, 
it appears safer to interpret its overall density profile as representing its dark matter 
distribution. This is probably because, the total stellar mass to the entire cluster 
mass is at most a few percent. So a galaxy in the cluster can raise this ratio locally, 
but the overall density profile is not significantly affected. 

4.2. Observational Status 

Enormous observational efforts have been made to obtain mass distributions 
of the central regions of galaxies and clusters of galaxies. Below we comment on 
the observational status from the point of view of our argument on the dynamical 
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stability. In many cases, the inner density power-law indices are poorly constrained 
due to finite spatial resolutions. As a practical division, we tentatively define a 
shallow-cusp/core profile for a < 1 as "a core" and a stable-cusp profile for 1 < 
a as "a cusp". This discrimination is based on our stability criterion discussed 
previously. We emphasize our previous caution that a flat-top stable core is often 
indistinguishable from a shallow unstable cusp. 

4.3. Centers of Galaxies 

The observational status is better off for spiral galaxies for which the abundance 
of interstellar medium allows measurements of rotation curves using Ha, CO, and 
HI 21cm emission lines. However, the limitation in spatial resolution often hampers 
conclusive division between a core and a cusp. 

For dwarf galaxies and low-surface-brightness galaxies (LSBs), the predominance 
of constant density cores have been claimed by many authors, 14 )~ 18 ) with varied 
confidence levels in their conclusions. For instance, Swaters et al. do not rule 
out cusp interpretations based on their analysis of data quality, 14 ) while Gentile et 
al. strongly claim the presence of a constant density core for one dwarf galaxy. 17 ) 
Nuclear activities or the presence of central black holes are generally not reported for 
dwarf galaxies and LSBs. If "a constant density core" implies "a flap-top profile", 
such a core is probably thermodynamically stable and a final product of dynamical 
evolution. 

Our Galaxy is a spiral galaxy and its rotation curve is well known. 19 ) At 0.1 < 
r < 20 kpc, the rotation curve is nearly flat, indicating an isothermal density profile 
of p(r) ~ r~ 2 . If our Galaxy were composed purely of collisionless dark matter 
particles, this isothermal profile could describe the entire Galaxy, since a = 2 implies 
dynamical stability. However in reality, the stellar mass fraction appears high in 
the inner galaxy, and dissipative collapse of baryons apparently happened in the 
past. The presence of a bulge and a black hole of 3 x 10 7 M Q at the Galactic center, 
indicates, star formation and subsequent black hole formation, near the center, a = 2 
is an isothermal profile and this may indicate the past relaxation of the mixture of 
baryons and dark matter on the scale of the entire Galaxy. 

The rotation curve should either rise (the presence of a black hole) or decline 
linearly (the presence of a constant core) towards the center, if it is measured with 
sufficient spatial resolution. Clearly both types exist for spiral galaxies in general. 20 ) 

For elliptical galaxies with little interstellar gas component, rotation curves are 
not available. However at least to our knowledge, the density profile has been ob- 
tained by the measurement of velocity dispersion of intracluster stars in the halo of 
a cD galaxy, A 2199. 21 ) They prefer a finite density core than a cusp for the dark 
halo density profile. For E/SO galaxies which happen to be the lenses of strongly 
lensed quasars, Treu and Koopmans have found that the total mass slope is about 
a = 1.75 ± 0.1, while the inner dark matter slope is a = 1.35, if the stellar ve- 
locity ellipsoid is isotropic, but a < 0.6, if the galaxies are radially anisotropic. 22 ) 
Observationally many ellipticals are expected to harbor black holes as dead quasars. 

In any case, the application of our stability criteria is limited for galaxies for the 
reasons mentioned above. 
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4.4. Centers of Clusters of Galaxies 

Gravitational lensing gives direct measurements of two-dimensional mass distri- 
butions of clusters of galaxies and thus is the most promising approach for obtaining 
volume density profiles under the assumption of spherical symmetry. Especially, 
strong lensing provides robust information on the two-dimensional encircled masses 
within tangential and radial critical radii, with minimal assumptions. Detailed 
strong lensing results obtained for a couple of nearby clusters, CL 0024+1654 23 ) 
and A1689, 24 ) indicate the presence of soft cores or at most shallow cusps, although 
there is a counter example of a steep dark matter profile in A383, which appears 
to be formed under the influence of its cD galaxy. 25 ' For A1689, Broadhurst et al. 
obtained the column density profile for the entire cluster by combining strong and 
weak lensing results. 26 ^ Their best-fit core-power-law profile indicates the presence of 
a finite density core. Sand, Treu and Ellis combined strong lensing and spectroscopy 
of gravitational arcs and the brightest cluster galaxy in the cluster, MS 2137-23, and 
analyzed the consistency of observations and the predicted cusps by N-body simula- 
tions. 27 ) They concluded that a < 0.9 with 99% confidence level and preferred a flat 
profile. The same group later obtained a mean power-law index a = 0.52, for three 
clusters with both radial and tangential arcs. 28 ) They suggest that the relationship 
between dark and visible matter in the cores of clusters is very complex. However, 
except for the case that the distribution of dark matter is heavily perturbed by that 
of luminous matter, it appears safe to consider that strong lensing results favor soft 
cores at the centers of clusters. 

The observational situations are summarized in Table ITTI 



Table II. Density profile observations. Here we define a core by a soft density profile with a < 1, 



and a cusp by a steeper profile with 


a > 1. 






Object Type 


Core/Cusp 


Method 


Comment 


Dwarf Galaxies 


Cores 


Rotation Curve 


No BHs 


Low Surface Brightness Galaxies 


Cores 


Rotation Curve 


No BHs 


Our Galaxy (Spiral) 


Cusp 


Rotation Curve 


BH 


Spiral Galaxies 


Cusps/Cores 


Rotation Curve 


With/ Without BHs 


E/SO Galaxies 


Cusps 


Lensing 


BHs as Dead QSOs 


Cluster Centers 


Cores 


Lensing 


No BH? 



§5. Concluding Remark 

In this paper, we have examined the requirements for the stability of spherically 
symmetric density profiles, first for ordinary gas with well-defined EOS, and then for 
collisionless particles, with a heuristic approach. A flat-top profile is thermodynam- 
ically stable near the origin, and appears to be a form of final product of evolution. 
In the case of a power-law density profile for a collision system, p(r) oc r~ a , we con- 
jecture that dynamical stability depends on a. The system is unstable and possibly 
in a transient stage toward a flat-top profile, for < a < 1. The system is thermo- 
dynamically metastable, but dynamically stable with a long life for 1 < a < 3. It is 
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unstable for a > 3, indicating a collapsing stage to the center. 

One physical way to support a flat-top density structure in a cluster of galax- 
ies, is to introduce degeneracy pressure of fermionic dark matter. In a compan- 
ion paper, 2 ) we model the density profile of the center of the cluster of galaxies, 
A1689, under the assumption that some fraction of dark matter is eV-mass degen- 
erate fermions/neutrinos. 
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